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@ Frequenzanalyseverfahren 

(g) Ein Frequenzanalyseverfahren umfaSt die Bewertung ei- 
nes zeitlichen, in Form von diskreten Abtastwerten vorlie- 
genden Eingangssignals mit einer Fensterfunktion. Anschlie- 
Qend wird das gefensterte Eingangssignal zur Erzeugung 
eines Satzes von Koeffizienten Fourier-transformiert. 
Um ein solches Verfahren so weiterzubilden, daS die 
Eigenschaften des menschlichen Ohres nicht nur hinsicht- 
lich der spektralen Abbildung in den Frequenzbereich, 
sondem auch hinsichtlich der Auflosung im Zeitbereich 
nachgebtldet werden, wird ein Block des Eingangssignals 
mit einem Satz von unterschiedlichen Fensterfunktionen zur 
Erzeugung eines Satzes von mit den jeweitigen Fensterfunk- 
tionen gewichteten SIdcken von Abtastwerten bewertet, 
deren Fourier-Transform ierte unterschiedliche Bandbreiten 
haben, bevor jeder der zeitgleich erzeugten BIdcke von 
Abtastwerten einer eigenen Fourier-Transformation derart 
" unterzogen wird, daft fur jede Fensterfunktion zumindest 
\ jeweils ein Koeffizient berechnet wird, dem die Bandbreite 
der Fourier-Transformierten dieser Fensterfunktion zugeord- 
net ist und da& die Koeffizienten so gewahit werden, da& die 
ihnen zugeordneten Frequenzbander im wesentlichen anein- 
ander angrenzen. 
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Beschreibung 

Die vorliegende Erfindung beirifft ein Frequenzanalyseverfahren, bei dem aus einem zeitdiskret. in Form von 
Blacken von Abtastwerten vortiegenden Eingangssignal mittels einer Fensterfunktion uberiappend aufeinander- 
5 foigende Bldcke von Abtastwerten entnommen werden, woraufhin die so mit den Fensterfunktionen bewerteten 
Bt6cke 2ur Erzeugung eines Satzes von Koeffizienten fourier-transformiert werden, nach dem Oberbegriff des 
Patentanspruchs 1. 

Es sind hauptsichlich zwei verschiedene Frequenzanalyseverfahren bekannt, mittels derer ein zeitliches, in 
Form von Abtastwerten vorliegendes Eingangssignal in eine Mehrzahl von Koeffizienten, d. h. in den Spektral- 

10 bereich transformiert werden kann. Diese Verfahren lassen sich in sogenannte Teilbandverfahren (vgL bspw. 
DE 40 1 1 577 Al) und Transformationsverfahren (vgL bspw. US 5t 42 656) kiassifizieren, wenngleich Transfor- 
mations- und Filterbankanalysen mathematisch zunSchst eng miteinander verwandt sind. 

FQr die Audiosignalanalyse und/oder -codierung ist es zweckmaBig und vorteilhaft, die Frequenzanalyse den 
Eigenschaften des menschlichen Ohres als Schallempfanger anzupassen. Eine wesentliche Erkenntnis der Psy- 

15 choakustik ist hierbei, daB das Ohr Schallereignisse, die in eine bestimmte Bandbreite, die sogenannte Frequenz- 
gruppenbreite fallen, gemeinsam bewertet Schallereignisse auBerhalb einer solchen Frequenzgruppenbreite 
werden getrennt wahrgenommen. 

Die Frequenzgruppenbreite ist — im Gegensatz zu rein technischen Frequenzanalysen — nicht konstant; mit 
zunehmender Mittenfrequenz wird auch ihre Bandbreite groBen Der Zusammenhang zwischen Frequenzgrup- 

20 penbreite und Mittenfrequenz ist daruberhinaus, wie in Fig. 2 gezeigt, nichtlinear. 

Dabei verwenden die sogenannten Teilbandverfahren zur Spektralanalyse Filterbanke, die eine spektrale 
Zerlegung entsprechend den Frequenzgruppen vornehmen. Oblicherweise verwendet man hierzu Kaskadenan- 
ordnungen aus Halbbandfiltem. Durch die fortgesetzte Halbierung des Basisbandes mittels einer derartigen 
Halbbandfilterkaskade kann nur eine grobe Annaherung an die Aufteilung in Frequenzgruppen erreicht werden. 

25 Die Zeitaufldsung von Bandfiltem ist umgekehrt proportional zu ihrer jeweiligen Bandbreite; sie wird also mit 
zunehmender Bandbreite feiner und ist damit gut den Eigenschaften des Gehors angepaBt. Halbband-Filterkas- 
kaden erzeugen aber, bedingt durch das Prinzip der fortgesetzen Halbierung des Basisbandes, relativ groBe 
Abweichungen der Frequenzaufteilung von der Vorgabe durch die Frequenzgruppenaufteilung, die Zwicker fiir 
das Gehdr gefunden hat. 

30 FFT-Analysen kdnnen zwar bei entsprechend feinem Frequenzraster eine genaue Aufteilung in Frequenz- 
gruppen gewahrleisien, die Zeitaufldsung ist dann jedoch zu grob und insbesondere konstant fur alle Frequenz- 
gruppen, was wiederum auch nicht den Eigenschaften des menschlichen Schallrezeptors entsprechen kann. wie 
die Wahrnehmbarkeit von Vorechost6rungen beweist 
Prinzipiell ist die Zeitaufldsung durch die Rekonstruktion einer Impulsfolge mit Hilfe eines idealen Tiefpasses 

35 gegeben. Die aus der Impulsantwort eines Filters ablesbare Zeitaufldsung wird als "resolution span" bezeichnet 
Sie ist in Fig. 1 wiedergegeben. 

Wie in Fig. 1 gezeigt ist, ist die Erzeugung der auf fg bandbegrenzten Funktion f(t) durch Superposition von 
mit dem jeweiligen Abtastwert bewerteten sinx/x-Funktionen an(t) moglich. Die sinx/x-Funktionen (Spaltfunk- 
lionen) stellen die Impulsantwort eines ideaien Tiefpasses mit der Grenzfrequenz fg dan Die Zeitaufldsung des 

40 idealen Tiefpasses ist T. 

Das Prinzip wird fiir ideale Interpolatoren verwendet Dazu werden sinx/x-Funktionen mit dem jeweiligen 
Abtastwert bewertet und um einen Zeitabschnitt T verschoben. der sich aus dem Abstand von Maximum zum 
ersten Nulldurchgang der sinx/x- Funktion ergibL Die Superposition aller dieser sinx/x-Funktionen ergibt die 
interpolierte Zeitfunktion. Der Veriauf der auf diese Weise rekonstruierten kontinuierlichen Zeitfunktion stimmt 

45 an den Nutldurchg^ngen der Spaltfunktionen mit den Abtastwerten exakt uberein. Diese Bedingung ist als 
Interpolattonsbedingung bekannt 

Als Zeitaufldsung wird eben dieser Abstand von Maximum zum ersten Nulldurchgang bzw. der Abstand der 
MulldurchgSnge der sinx/x-Funktion bezeichnet Es ist leicht einzusehen, daB die Zeitaufldsung gleich der 
Abtastperiode ist, wenn als Interpolationsfilter ein TiefpaB mit einer Grenzfrequenz fg gleich der halben 

50 Abtastfrequenz fa. also ein AllpaB gewahit wird. Das Produkt aus Bandbreite und Zeitaufldsung betrSgt in 
diesem Fall B-T » fa - */2fa » 0,5. Bei einer bestimmten Bandbreite kann die damit gegebene Zeitaufldsung nicht 
unterschritten werden. Dies wird in Analogic zu einem ahnlichen, aus der Physik bekannten Sachverhalt als 
Heisenberg-Grenze bezeichnet Teilband- und/oder Transformationsverfahren zur spektraien Zerlegung fmden 
sich in fast alien heute bekannten Codern. Beide sind fiir eine gehdrgerechte ReprSsentation im Frequenz- und 

55 Zeitbereich ungeeignet Speziell muB ein MeBverfahren, das beansprucht, die Qualitat heutiger Codierverfahren 
beurteilen zu kdnnen, sicher bereits hinsichtlich der Frequenzanalyse genauer arbeiten; erst dann ist die Genau- 
igkeit, mit der hieruber die Mithdrschwellen bestimmt werden hdher, als in den zu beurteilenden Coder-Deco- 
der-Kombinationen. 

Die Nachteile herkdmmlicher, heute fast ausschlieBlich angewendeter Verfahren zur gehdrrichtigen Fre- 
60 quenzanalyse seien in zusammengefaBter Form verdeutlicht: 

— Teilbandverfahren arbeiten zwar an der Heisenberg-Grenze; die Frequenzgruppenbreiten werden 
jedoch nur grob getroffen. 

— Transformationsverfahren sind zwar in der Lage, die Frequenzgruppenbreiten im Prinzip beliebig genau 
65 nachzubilden, die Zeitaufldsung ist jedoch bei weitem zu grob, dariiber hinaus konstant und somit nicht 

gehdrad^quat 

— Hybride Frequenzanalysen nahem sich der geforderten Zeit-/Frequenzaufldsung des Gehdrs zwar 
besser. eine exakte Nachbildung ist aber auch hiermit prinzipiell nicht moglich. 
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Der Erfindung liegt die Aufgabe zugrunde, ein Frequenzanlayseverfahren zu schaffen, das die Eigenschaften 
des menschlichen Ohres nicht nur im Frequenzsondem gleichzeitig auch im Zeitbereich nachbildet Diese 
Aufgabe wird durch ein Frequenzanalyseverfahren geraaO Anspruch 1 gelosL 

GemaB der Erfindung werden bei einem Frequenzanalyseverfahren der eingangs genannten Art aus einem 
zeitlich in Form von Abtastwerten vorliegenden Eingangssignai mitteis Fensterfunktionen Qberlappend aufein- 5 
anderfolgende BIdcke von Abtastwerten entnommen, worauf die so mit den Fensterfunktionen bewerteten 
Blocke zur Erzeugung eines Satzes von Koeffizienten Fourier-transformiert werden, wobei aber jeweils ein 
solcher Block des Engangssignales mit einem Satz von unterschiedlichen, symraetrischen Fensterfunktionen 
gieicher Lange zur Erzeugung eines Satzes von mit den jeweiligen Fensterfunktionen gewichteten BIdcken von 
Abtastwerten bewertet wird Die Fourier-Transformierten der Fensterfunktionen besitzen, bedingt durch die to 
unterschiedlichen Fensterfunktionen, auch unterschiediiche Bandbreiten. Zu jedem der zeitgleich erzeugten 
BIdcke von gefensterten Abtastwerten liBt sich uber eine eigene Fourier-Transformation mindestens ein Koeffi- 
zient berechnen, dem die Bandbreite der Fourier-Transformierten dieser Fensterfunktion zugeordnet ist, wobei 
die Koeffizienten so gewahJt werden, daQ die ihnen zugeordneten Frequenzb§nder, bezogen auf ihre Bandbreite, 
sich nur geringfflgig Qberiappen oder beabstandet sind. Mit anderen Worten schlieBen die den Koeffizienten 15 
zugeordneten Frequenzbander im wesentiichen luckenlos und uberdeckungsfrei aneinander an. 

Ein bevorzugtes Ausfuhrungsbeispiel des erfindungsgemaBen Frequenzanalyseverfahrens wird nachfolgend 
unter Bezugnahme auf die beiliegenden Zeichnungen naher erliutert Es zeigen: 

Fig. I eine Darstellung zur Verdeutlichung. daB die Erzeugung einer auf Fg bandbegrenzten Funktion f(t) 
durch Oberlagerung von mit dem jeweiligen Abtastwert bewerteten sinx/x- Funktionen mdglich ist; 20 

Fig. 2 den Zusanmienhang zwischen der Frequenzgruppe AFg in Abhangigkeit von der Mittenfrequenz fm; 

Fig. 3 den Zusammenhang zwischen einer Frequenzachse in Hertz ais Ordinate und einer Frequenzachse in 
BARK ais Abszisse; 

Fig. 4 und 5 Frequenz- und Zeitdarsteilungen zur Verdeutlichung des Ahnlichkeitssatzes der Fourier-Trans- 
formation zur Verdeutlichung, daB zu einer Verbreiterung der Bandbreite eines Spektralkoeffizienten eine 25 
Verschmalerung der Fensterfunktion im Zeitbereich natig ist; 

Fig. 6 ein Beispiel fur eine Fourier-Transformierte einer Fensterfunktion, welche die Energieerhaltung ge- 
wahrleistet; 

Fig. 7 ein Beispiel fur eine Fourier-Transformierte, die eine feinere Zeitaufidsung. jedoch nur eine reduzierte 
Frequenzaufldsungbesitzt; 3q 

Fig. 8 ein FluBdiagramm eines BARK-Generators zur Ermittlung von Fensterfunktionen im Zeitbereich, die 
frequenzgruppenbreite Spektralkoeffizienten erzeugen; 

Fig. 9 eine Darstellung zur Beurteilung eines BARK-Generatorentwurf es; 

Fig. 10 eine Darstellung des iterattonsergebnisses des BARK-Generators in Abhangigkeit von der Rasterwei- 
te; - 35 

Fig. 1 1 eine Darstellung des Quotienten aus der Anzahl unterschiedlicher Frequenzgruppen und der Gesamt- 
zahl der Bander in Abhangigkeit von der Rasterzahi; 

Fig. 12 ein Blockdiagramm einer Vorrichtung zur Durchfuhrung der BARK-Transformation auf der Grundla- 
ge einer 512-Punkte-FFr; 

Fig. 1 3 ein rekursives Filter erster Ordnung zur Bildung des Goertzel-AIgorithmus; 40 

Fig. 14 ein rekursives Filter zweiter Ordnung zur Bildung eines Goertzel-AIgorithmus zweiten Grades; 

Fig. 15 eine transponierte Struktur eines Filters zur Bildung eines Goertzel-AIgorithmus zweiten Grades; 

Fig. 16 eine bevorzugte Ausfuhrungsform eines rekursiven Filters zur Bildung des Goertzel-AIgorithmus 
zweiten Grades ais Vorrichtung zur Berechnung eines Koeffizienten der BARK-Transformation; 

Fig. 17 eine Darstellung zum Vergleich des Rechenaufwandes fur einen komplexen Spektralkoeffizienten bei 45 
Verwendung einer schnellen Fourier-Transformation gegenuber dem Goertzel-AIgorithmus. aufgetragen uber 
verschiedene Transformationslangen N; 

Fig. 18 eine Darstellung der Anzahl GA der Goertzel-Algorithmen, die innerhalb des Rechenzeitbedarfs fQr 
eine schnell Fourier-Transformation auf einem Signalprozessor gerechnet werden konnen, aufgetragen uber 
verschiedene Transformationslangen N; 5q 

Fig. 19 eine Darstellung der Betragsubertragungsfunktion der Einzelfilter einer BARK-Transformation; und 

Fig. 20 eine Darstellung des Frequenzganges ais Summe der Amplitudenubertragungsfunktionen. 

Transformationsverfahren lassen sich offensichtlich mit deutlich weniger Aufwand ais Filterbanke realisieren. 
Auf der anderen Seite sind sie aber nur dann mit Erfolg einsetzbar. wenn auf der spektralen Achse (der 
Frequenzachse in Hertz) lediglich aquidistant angeordnete Reprasentanten benCtigt werden. Eine nicht aquidist- 55 
ante Anordnung auf der spektralen Achse uber eine Gruppierung herzustellen. fuhrt zu einer Realisiening, die 
von vome herein nicht gehorrichtig ist. weil sie nicht entlang der Heisenberg-Grenze arbeitet 

Der Zusammenhang zwischen einer Frequenzachse in Hertz und einer Frequenzachse in BARK ist aber — 
jedenfalls fQr niedrige Frequenzen — naherungsweise linear. Fig. 2 [Zwicker, E; Feldtkeller, R.: Das Ohr ais 
Nachrichtenempfanger, Z AufU Hirzel Veriag, Stuttgart 1967] laBt dies bereits erkennen. Deutlicher wird der eo 
lineare Zusammenhang, wenn man die Frequenz in Hertz Uber der Tonheit in BARK auftragt, wie Zwicker dies 
in Fig. 3 getan hat [Zwicker, E; Feldtkeller, R.: Das Ohr ais Nachrichtenempfanger, 2. Aufl., Hirzel Veriag. 
Suttgart 1967]. Dann ist erkennbar, daB der Zusammenhang zwischen einer Frequenzachse f in Hertz (Ordinate) 
und einer Frequenzachse z in BARK (Abszisse) zumindest im unteren Frequenzbereich bis etwa 5 BARK linear 
ist; das bedeutet, das hier die Teilbander gleich breit sind Im Bereich von 0 Hz bis etwa 500 Hz kann also in guter 55 
Naherung von einem linearen Zusammenhang gesprochen werdea In diesem Bereich laBt sich daher fur jede 
Mittenfrequenz in etwa die gleiche Frequenzgruppenbreite von 100 Hz angebea Hier wird zunachst eine 
Fourier-Transformation vom Zeitbereich in den Frequenzbereich eingesetzt Sie kann hier ohne die Qbliche 
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Gruppierung von Spektralwerterv d. h. mit einem idealen Produkt von Bandbreite und Zeitauf!6sung verwendet 
werden. Man gewinnt dabei den Vorteil, zumindest die ersten funf Spektralkoefflzienten gemeinsam in einem 
Ablauf (z. B. in einer schnellen Diskreten Fourier-Transformation) berechnen zu kannen. Oblicherweise ergeben 
sich bei der Abbildung einer Frequenzachse von 0 Hz bis etwa 22 kHz auf eine BARK-Achse etwa 25 Spektralli- 

5 nien. Sie kdnnen als Repr^sentanten fur die frequenzgruppenbreiten Binder angesehen werden. Die Berech- 
nung der ersten fQnf Spektrallinien mittels einer Fourier-Transformation hat mit einer Schrittweite auf der 
Frequenzachse zu erfolgen, die der in diesem Bereich ann^emd konstanten Frequenzgruppenbreite, also etwa 
too Hz, entspricht Bei einer Abtastfrequenz von beispielsweise 48 kHz muO demnach mit einer 480-Punkte 
Fourier-Transformation gerechnet werden. Von den 240 Spektrallinien kdnnen dann schon die ersten 5 verwen- 

iO det werden, die restlichen 235 Spektrallinien mflssen verworfen werden. da sie zu hSheren Frequenzen hin keine 
Frequenzgruppenbreite mehr reprSsentieren (ihre Bandbreiten sind dafur zu schmal). FQr den geforderten 
Entwurf entlang der Heisenberg-Grenze sind wegen der Reziprozitgt von Bandbreite und Zeitaufl6sung B-T 
0,5 zur Erzielung bretterer Binder entsprechend schmalere Fenster im Zeitbereich fOr weitere, getrennt zu 
berechnende Fourier-Transformationen zu wihlen; diese ergeben dann — wie gewflnscht — auch eine bessere 

15 Zeitaufldsung. 

Der Ahnlichkeitssatz der Fourier-Transformation besagt, daB zu einer "Dehnung" der Bandbreite eines 
Spektralkoeffizienten eine "Stauchung"der Fensterfunktion im Zeitbereich ndtig isL Zeitaufldsung T und Band- 
breite B verhalten sich in reziproker Weise zueinander. Diese Reziprozitat wird mittels der in Fig. 4 und Fig. 5 
veranschaulichten Korrespondenzen der Fourier-Transformation gezeigt Zur Erzielung einer hdheren Band- 

20 breite muB demnach die Fensterfunktion im Zeitbereich gestaucht werden; seine Zeitaufldsung wird damit 
besser. Der dem erfindungsgemaOen Frequenzanalyseverfahren zugrunde liegende Gedanke ist es also, die 
Ordnung der Fourier-Transformation (die Anzahl der Koeffizienten) beizubehaiten, gleichzeitig aber dafQr zu 
sorgen, daB die hdheren Spektrallinien Reprisentanten fur breitere Bander darstellen. Auf diese Weise erreicht 
man auch gleiche Steiiheiten fQr die Filterflanken. Dies ist die Voraussetzung fOr ein konstantes GesamtQbertra- 

25 gungsverhalten fflr den Betragsfrequenzgang, das fur den Einsatz in MeBverfahren, insbesondere aber auch in 
Codierverfahren, gefordert werden muB. 

Mit einer BARK-Transformation wird eine gehorgerechte Transformation vom Zeitbereich in den Frequenz- 
bereich bezeichnet, wobei eine aquidistante Darsteilung der Spektralkoeffizienten auf einer BARK-Achse (und 
nicht auf einer Hertz- Achse) erreicht wird. Die hier entwickelte BARK-Transformation unterscheidet sich von 

30 anderen (meist nur gehdrangepaBten) Abbildungen auf die BARK-Achse dadurch, daO sie entlang der Heisen- 
berg-Grenze und damit gehdrrichtig arbeitet 

Eine Fourier-Transformation, die grdBere -6 dB- Bandbreiten liefert, als es dem Abstand Af der Spektralli- 
nien entspricht, ist im herkdmmlichen Sinne naturlich unbrauchbar (Fig. 6). 
Normalerweise ist die Bandbreite Bn aller Spektralkoeffizienten der Fourier-Transformation durch die Ord- 

35 nung N und die Abtastfrequenz fa vorgeschrieben: Bn =» fa/N. Die Betragsubertragungsfunktionen schneiden 
sich dann in |H(f)| » 0,5, d. h. an den -6 dB-Punkten (Fig. 7). Dennoch kdnnen von einer nach Fig. 6 modifizier- 
tcn Fourier-Transformation nur die Linien fur die BARK-Transformation verwendet werden, deren Bandbreiten 
bzw. Mittenfrequenzen in das von Z wicker fiir die Frequenzgruppenbreiten gegebene Schema passen. Durch die 
— zumindest innerhalb eines Rasters — freie Wahl der Bandbreite jedes Koeffizienten erubrigt sich eine 

40 Gruppierung von einigen Koeffizienten zu einer Frequenzgruppe mit all ihren Nachteilen, wie sie bei herkdmm- 
lichen gehdrangepaBten Transformationsverfahren ndtig ist Weil das Frequenzraster fur jede neu zu suchende 
Linie einer BARK-Transformation beibehalten wird, lassen sich nun BSnder aneinanderreihen, deren Breite in 
beliebig guter NSherung der angestrebten Abbildung auf die BARK-Achse entspricht Ein Entwurf entlang der 
Heisenberg-Grenze ist dabei gewahrieistet Die Staffelung ist deshalb luckenlos" mdglich, weil Bandbreiten 

4S gewahlt werden, die genau in das durch die LSnge der Fourier-Transformation gegebene Af-Raster passen. Die 
Summe aller AmplitudenQbertragungsf unktionen bieibt dabei dennoch konstant: die Flankensteilheit der bei der 
Staffelung entstehenden Bandfilter ist nur von der Anzahl der Filterkoeffizienten (der FensteriSnge), nicht aber 
von der Bandbreite abhangig. Diese Anzahl ist aber konstant, sie entspricht normalerweise (d. h. ohne Polypha- 
senaddition) der LSnge der Fourier-Transformation. Weil alle Flanken die gleiche Steilheit aufweisen, treten 

50 unerwQnschte Schwankungen im Gesamtiibertragungsverhalten (Summenf requenzgang) hier nicht auf. 

Die in Form von Abtastwerten vorliegenden Fensterfunktionen mussen nicht notwendigerweise die gleiche 
Anzahl von Abtastwerten aufweisen, die der LSnge der auf die Fensterung folgenden Fourier-Transformation 
entspricht Vielmehr kann durch eine Veriangerung der Fensterfunktion uber die Transformationsiange hinaus 
eine Verbesserung der Flankensteilheit der Einzelfilter erreicht werden. Das Additionstheorem der Fourier- 

55 Transformation ertaubt es, auch llngere Sequenzen als durch die Transformationsiange gegeben in den Fre- 
quenzbereich zu transformieren. Hierzu wird eine Anzahl von aufeinanderfolgenden Teilbldcken aus dem 
bereits gefensterten langeren Engangsblock von Abtastwerten gebildet nach dem Prinzip der Superposition 
addiert und gemeinsam in den Frequenzbereich transformiert Dieses Prinzip ist auch als Polyphasen-Filterung 
bekannt[Chrochiere, R.E.; Rabiner L: Multirate Digital Signal Processing, Prentice-Hall, Engtewood Cliffs, New 

60 Jersey, 1983]. Durch die Polyphasenfilterung wird insbesondere die Zeitaufidsung nicht bertlhrt, da die dem 
Filterentwurf zugrunde liegende Spaltfunktion die gleiche bieibt 

Die Erfindung liefert femer ein als BARK-Generator bezeichnetes Iterationsverfahren, das nunmehr erlautert 
wird. Die Staffelung der Bander wird nach Vorgabe einiger Parameter einem Rechenautomaten Qberiassen. 
Diese Staffelung liefert fUr jede Frequenzgruppe einen Satz von Filterkoeffizienten fiir das Zeitfenster einer 

65 schnellen Fourier-Transformation. Es hat sich als gOnstig erwiesen, die Staffelung ab der halben Abtastfrequenz 
(d. h. der hdchsten vorkommenden Signalfrequenz) zu tieferen Frequenzen hin und nicht umgekehrt vorzuneh- 
men. (Eine Staffelung ab der niedrigsten im Signal vorkommenden Frequenz ist natQrlich auch mdglich). Die 
Staffelung gewahrleistet nicht von vome herein, daB am Ende die untere Grenzfrequenz der untersten Fre- 
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quenzgruppe mit 0 Hz zusammenfilit Brauchbar sind aber auch f terationen. die untere Grenzfrequenzen fur das 
unterste Band bis hinauf zu etwa 50 Hz Uefem. Ein mdglicherweise im Signal enthattener. normaJerweise 
unhorbarer Gleichanteil wird dadurch in die Messung stark gedampft eingehen. Audtosignale reichen nicht bis 
zu einer Frequenz bis 0 Hertz hinab. Im Verlauf der Iteration werden also wettere B^der- unter die bereits 
ermittelten gestaffelt Die Bandgrenzen der frequenzgnippenbreiten Bander, aber auch deren Mittenfrequenzen 5 
fallen dabei auf ein vorzugebendes Raster. 

Bei der Suche nach den Koeffizientensatzen fur weitere Fensterfunktionen wird nach dem in Fig. 8 dargestell- 
ten RuOdtagramm vorgegangen{CBW = Frequenzgruppenbreite nach Zwicker). 

— Die obere Bandgrenze ist bekannt Sie ergibt sich aus der unteren Bandgrenze des bereits ais passend 10 
ermittelten Zeitfensters fur die vorherige, hdhere Frequenzgruppe, wobei 6 dB-Bandbretten zugrunde 
gelegt werden; bei Iteradonsbeginn wird die obere Grenzfrequenz auf die halbe Abtastfrequenz als Start- 
wert gesetzt 

— Die Mittenfrequenz des als nichstes zu suchenden Bandes wird nun solange auf dem Raster verschoben, 
das sich durch die Abtastfrequenz und die Ordnung der Fourier-Transformation ergibt^ bis aus dem 15 
Ensemble von moglichen Bandbreiten eine gefunden wird, die der Frequenzgruppenbreite von Zwicker am 
n^chsten kommt Ist das verwendete Frequenzraster nur fein genug, so wird auch die Genauigkeit, mit der 
eine Frequenzgruppenbreite getroffen wird, beliebig gut sein. 

— Aus der gefundenen Bandbreite ergibt sich der Koeffizientensatz fur eine der Fensterfunktionen und 
damit die Zeitauflosung dieses Fensters gemaB der Heisenberg-Grenze, d h. die Anfordening im Frequenz- 20 
bereich bestimmt dabei den Verlauf der Fensterfunktion im Zeitbereich* Hier wird das Filterentwurfsver- 
fahren nach Tietze, U.; Schenk. Ch-: Halbleiterschaltungstechnik, 9. Auflage, Springer- Verlag, Berlin, 1989, 
eingesetzt 

— Die untere Grenzfrequenz des gefundenen Bandes bildet die obere Grenzfrequenz des als n^chstes zu 
ermittelnden Bandes. Wird im Verlauf der Iteration die Mittenfrequenz negativ, so ist das Abbruchkriterium 25 
erreicht 

Es ist uberraschend, daB bei der Untereinanderstaffelung der frequenzgruppenbreiten Binder weit weniger 
als 25 unterschiedliche Zeitfenster notig sind Dies ist einerseits auf die nahezu konstante Frequenzgruppenbrei- 
te zu tiefen Frequenzen bis etwa 5 BARK zurilckzufuhren (vgl. Fig. 3); andererseits ist aber auch entscheidend, 30 
daB durch die Vorgabe eines Frequenzrasters eben nur diskrete Mittenfrequenzen bzw. Bandgrenzen eriaubt 
sind Wie Fig. 9 zeigt, konnen bereits bei einer Diskretisierung auf 512 Frequenzen bis zur halben Abtastfre- 
quenz die 25 frequenzgruppenbreiten Bander ausreichend genau mit nur 14 verschiedenen Fensterfunktionen 
nachgebildet werden: Die Beurteilung des BARK-Generator-Entwurfs ist sofort graphisch am Bildschirm mog- 
lich. Die Abweichung der Bandbreiten von den Frequenzgruppenbreiten nach Zwicker ist bereits bei einem 35 
reiativ groben Frequenzraster von 512 diskreten Frequenzen ausreichend klein (Beispiel in Fig. 9 fur 44,1 kHz 
Abtastfrequenz). 

Der BARK-Generator leistet die foigenden Aufgaben: 

— Der Entwurf ist prinzipiell fiir jede Abtastfrequenz moglich, speziell aber fiir die gangigen Abtastfre- 40 
quenzen 44,1 kHz, 48 kHz und 32 kHz; eine Abtastfrequenz von 44.1 kHz ist fur Compact- Disc Musikmate- 

rial ublich; Abtastfrequenzen von 32 kHz und 48 kHz findet man bei DAT-Medien (Digital Audio Tape). 

— Untereinanderstaffelung von frequenzgruppenbreiten Bandern. Die Genauigkeit, mit der der Entwurf 
den Frequenzgruppenbreiten von Zwicker entspricht, laBt sich durch die Vorgabe eines Rasters steuern. 

— Zu Versuchszwecken lassen sich auch Bandbreiten fur die Teilftlter einstellen, die Bruchteilen von einem 45 
BARK entsprechen: dies erdffnet die Mdglichkeit, das Tonh6henaufl6sungsvermdgen des menschlichen 
Ohres von etwa 1/27 BARK nachzubilden. 

~ Berechnung und Ausgabe der Satze von Fiiterkoeffizienten fiir die bendtigten Fenster unter Bealcksich- 
tigung der Polyphasenaddition in verschiedenen Ausgabeformaten. Eine der Ausgabemdglichkeiten ist eine 
Konstantentabelle, die sofort als Quelltext in ein Programm zur BARK-Transformation eingebunden 50 
werden kann. Die Bezeichnung BARK-Generator ist damit gerechtfertigt 

~ Bandbreite und Zeitauflosung der Fensterfunktionen ergeben (^nlich wie bei der trivialen Bandfilter- 
Kaskade) immer einen Entwurf entlang der Heisenberg-Grenze, 

— Eine visuelle Beurteilung der Iterationsergebnisse ist sofort durch die eingebundene Graphikausgabe 
mdglich. Eventuelle Modifikationen konnen interaktiv erfolgea 55 

Wie Fig. 10 zeigt l^uft die Iteration nicht fur beliebige Frequenzraster, sondem nur bei einigen ausgezeichne- 
ten Rasterweiten erfolgreich ab. Die Iteration kann immer dann als erfolgreich gelten, wenn die untere Grenz- 
frequenz des untersten Bandes im Bereich von etwa 0 bis 50 Hz liegt Es ist giinstig, einen im Signal mdglicher- 
weise enthaltenen Gleichanteil, der Ja normalerweise unhorbar ist, auch in die Messung stark gedimpft eingehen eo 
zu lassen. Der DurchlaBbereich der untersten Frequenzgruppe muB daher nicht unbedingt bis 0 Hz hinabrei- 
chen. Die Gesamtzahl der Bander, in die der Frequenzbereich von z. B.O bis 22.05 kHz unterteilt wird nahert sich 
mit zunehmender Rasterzahl (d h- mit feiner werdendem Frequenzraster) der Frequenzgruppenzahl von 25. 
Diese Anzahl ist fur die in diesem Beispiel gewShlte Abtastfrequenz richtig. 

Fig. 10 zeigt auch, daB die Anzahl unterschiedlicher Fensterfunktionen dabei inuner unterhalb der gesamten es 
Anzahl von Frequenzgruppen bleibL Hierin liegt ein wesentlicher Vorteil des BARK-Generators: Es sind — fur 
alle praktisch relevanten Fille — weniger Bewertungen des Eingangsvektors des Zeitsignals mit einer Fenster- 
funktion nOtig als durch die Transformation frequenzgruppenbreite Bander erzeugt werden. Bildet man den 

5 



BNS0OaD:<0E 4316297C1> 



DE 43 16 297 CI 



Quotienten aus der Anzahl unterschiedlicher Frequenzgnippen und der-Anzahl der Binder, so ergibt sich fur 
eine Rasterzahl von etwa 512 ein Minimum (Fig. 1 1). 

Hier arbeitet die BARK-Transformation am okonomischsten: Im Verbal tnis wird die geringste Anzahl unter- 
schiedlicher Binder bendtigt Noch kleinere Rasterzahlen sind zwar mdglich, fuhren in diesem Beispiel aber zu 
5 Abweichungen von der Frequenzgnippen breite nach Zwicker, die nicht mehr befriedigen. 

MaBnahmen zur Reduktion des Rechenaufwandes 

Der Rechenaufwand zur Ermittiung einer vollstandigen BARK-Transformation ist noch relativ grofl, da fOr 
10 jedes unterschiedliche Fenster eine neue Diskrete- Fourier-Transformation gerechnet werden mu6. Diese Vor- 
gehensweise geht mit einem meist nur begrenzt zur VerfQgiuig stehenden Potential an Rechenletstung offen- 
sichtlich sehr undkonomisch um, da die meisten der berechneten Fourier- Koeffizienten verworfen werden; nur 
sehr wenige (Coeffizienten (in vielen F^len sogar nur ein einziger) werden ffir die BARK-Transformation 
weiterverwendet 

IS Wie in Fig. 12 angegeben, kann eine BARK-Transformation, deren Grundlage eine 512-Punkte FFT bildet, 
bereits mit nur 1 1 unterschiedlichen Teilfenstem berechnet werden. Bei der Berechnung der Koeffizienten fflr 
die 25 frequenzgruppenbreiten Teilfilter uber eine schneile Fourier-Transformation (FFT) mflssen aber 99^6% 
aller FFT-Koeffizienten verworfen werden, nur 25 aus 5632 berechneten Koeffizienten werden tatsachlich 
verwendet. 

20 Fig* 12 verdeutlicht die Verschwendung von RechenkapazitSt in einem Beispiel mit 11 unterschiedlichen 
Fenstem und einer Rasterung auf 512 diskrete Frequenzen bis zur Abtastfrequenz: FQr die untersten 10 
Spektrallinien kann wegen des in diesem Frequenzbereich nSherungsweise linearen Zusammenhangs zwischen 
Frequenzgruppenbreite und Mittenfrequenz noch ein und dasseibe Zeitfenster verwendet werden. Fur die 
folgenden 3 Binder (U bis 13) wird ein neues Zeitfenster und damit eine neue gemeinsame FFT benfitigt Dies 

25 gilt auch noch filr die darauf folgenden 2 Frequenzgnippen; erst ab der I6ten Frequenzgruppe aufwSrts muB 
jeweils ein neues Zeitfenster pro Frequenzgruppe eingesetzt werden. Bei der ersten FFT werden dadurch 502 
von 512 Koeffizienten verworfen, dann 509 von 512. 510 von 512 und schiieBlich 8 mal 51 1 von 5\Z Insgesamt 
werden von 1 1»512 = 5632 Koeffizienten nur 25 verwendet Es mussen also 5607/25 ^ 224 mal mehr komplexe 
Fourier- Koeffizienten berechnet werden, als dann fflr die BARK-Transformation gebraucht werden. Wie in 

30 Fig, 12 jedenfalls beispielhaft skizziert wird, ergabe die Berechnung der BARK-Transformation iiber eine 
Anzahl von schnellen Fourier-Transformationen mit jeweils unterschiedlichen Zeitfenstem einen unverhaitnis- 
m^Qig hohen Rechenaufwand, wenn nicht eine efflzientere Methode zur Berechnung einzelner Koeffizienten als 
die schneile Fourier-Transformation eingesetzt wurde. 
Eine triviale M6giichkeit, nur einen Koeffizienten einer Fourier-Transformation zu berechnen, besteht darin, 

35 die Bestimmungsgleichung f Qr eine diskrete Fourier-Transformation anzuwenden. Mehr Vorteile bietet aber der 
Einsatz des Goertzel-AIgorithmus zu Berechnung einzelner Linien. Mehr Vorteile bietet aber der Einsatz des 
Goertzel-AIgorithmus zur Berechnung einzelneriinien [Goertzel, G.: An Algorithm for the Evaluation of Finite 
Trigonometric Series. Amer. Math. Monthly 65 (1 968), S34— 35]. 
Ausgehend von der z-Transformation X(z) einer Sequenz x(n) von Abtastwerten der L^nge N 

40 

X(z)=xW)-2» 

45 

ergibt sich die Diskrete Fourier-Transformation (DFT) mit 

Z « « (e-i2n/N)k „ e-j2Kk/N 

50 Ober k werden diskrete Frequenzen ausgewihlt 

Mit der Substitution z" = W^*^ = g-j2nk7N ^^^^^ Qpj ^m^j^ Polynom angegeben werden. Beispielsweise 
ergibt sich ffir N « 4 ein Polynom dritter Ordnung 

X(z) « x(3).z3 + x(2).z2 -h x(l).z x(O). 

55 

Ober das Horner-Schema 16st man die Multiplikationen mit Potenzen von z auf 

X(z) - ((x(3)-z + x(2)).z + x(I)).z -h x(0). 

60 Unter Einfilhrung einer neuen Variablen y(n) kann daraus eine lineare Differenzengleichung angegeben 
werden 

y(n) = z-y(n-l) + x(N-nX 

65 Anfangsbedingung ist hierbei y(0) 0. Das gesuchte Ergebnis fiir die Fourier-Transformierte findet man fOr n 
« N 

X(z) = y(N). 

6 • 
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Die lineare Differenzengleichung kann auch als rekursives Filter erster Ordnung aufgefaBt werden. Die 
eingangs gemachte Substitution z" = W"'' wird hierzu wieder ruckgangig gemacht 

y(n) - W*'.y(n-l) + x(N-n);y(0) = 0. 

5 

Dte damit gegebene Anordnung ist in Fig. 13 wiedergegeben. 

Die Eingangssequenz muO allerdings in umgekehrter Reihenfolge vortiegea 

AuQerdem wirkt sich nachteilig auf die Stabilttat and den Rechenaufwand fQr das Filter aus, daB die Zustands- 
variable y(n-l) fiber einen komplexen KoefRzienten W*^ auf den Filtereingang zuruckgefuhrt wird. 

Die Anzahl der Multiplikationen und Additionen ist zwar die gleiche wie bei der direkten Berechnung der lo 
diskreten Fourier-Transformation, die Drehf aktoren mussen bier jedoch nichi explizit berechnet oder fiber eine 
Drehfaktortabelle ("twiddle table") zur Verffigung gestellt werden; durch die rekursive Berechnungsweise 
werden sie vielmehr in jedem Durchlauf, ausgehend von einem Anfangswert, fur einen neu hinzukommenden 
Abtastwert der Engangssequenz x(n) weitergedreht 

Die Anzahl der bendtigten Multiplikationen laBt sich auf die Halfte reduzieren, wenn man die Euler'sche 15 
Gleichung einsetzt 

W*^ « e-j2«/N « cos(27ik/N) - j sin(2nk/N) 

und auf eine indizierte Notation fibergeht 20 

yn " W'Vn-l + XN-n =» [C0S(2nk/N) - j sin(27tk/N)]-yn.| + XN-m 

wobeiyn-i « [cos(2;ik/M) — j sin(27ik/N)] yn-2 + XN-n+t- 
Nach einigen weiteren Umformungen ergibt sich schlieBlich fur yn 25 

yn - [2 • cos(2^k/N)] • yn-i + yn-2 = -XN-n + [cos(2^k/N) + jsin(2|ik/N)] • XN+i-n. 

Durch Vergleich mit der allgemeinen Differenzengleichung zweiter Ordnung 

30 

y(k+2) + c, •y(k + l) + co-y(k)-b2-x(k + 2) + b, • x(k+ 1) + bo • v(k) 

kdnnen die Filterkoeffizienten fur ein rekursives Filter zweiten Grades in der ersten kanonischen Grundform 
angegeben werden, das sich wie in Fig. 1 4 gezeigt, darstellen laBt 

35 

ci = — 2cos(2nk/N);co = 1 

b2 - 0; bi = (27ik/N) + j sin (2nk/N); bo = 1. 

In dieser Form besttzt die Gleichung reelle FCoeffizienten fur die RQckfuhrung, allerdings einen komplexen 
Eingangskoeffizienten. 40 

Das rekursive Filter zweiten Grades in der zweiten kanonischen Struktur kann — wie in Fig. 15 gezeigt — 
noch transponiert werden. 

Da die Etngangsfolge rein reell ist, aber in umgekehrter Reihenfolge am Eingang anliegt, gilt ffir den ungera- 
den Anteil der Eingangsfolge 

hu(t)=-h(t). 

Der ungerade Anteil hu(t) liefert aber nur den ungeraden Imaginarteil der Fourier-Transformation; speist man 
die Eingangsfolge in nicht umgekehrter Reihenfolge ein, muB der Imaginarteil der Fourier-Transformierten 
lediglich negiert werden. 50 

Fig, 16 zeigt die endgultige Version des Goertzel-Algorithmus. Die Eingangsfolge x(n) liegt jetzt in normaler 
Reihenfolge am Eingang an. Der Ausgangswert muB erst bestimmt wenn der Eingangsvektor x(n) vollstandig 
abgearbeitet ist Zur Berechnung einer Spektrallinie ist es vorteilhaft zunachst aile Werte der Eingangsfolge x(n) 
in das Filter "hineinlaufen" zu lassen und erst dann den Ausgangswert y(n) fiber eine komplexe Multiplikation aus 
den Zustandsvariablen Ai und A2 zu berechnen. 55 

Mit diesem wenig bekannten Algorithmus ist es also moglich, einzelne Linien einer diskreten Fourier-Trans- 
formation zu berechnen; die dabei verwendete Struktur laBt sich auch als rekursives Filter auffassen. Ein 
besonderer Vorteil gegenfiber der FFT ist, daB hier keine Drehfaktortabelle bendtigt wird. AuBerdem treten 
innerhalb des Filters lediglich rein reelle Additionen und Multiplikationen auf. Nur am Filterausgang wird eine 
komplexe Multiplikation benotigt; zur Berechnung eines BARK-Koeffizienten muB diese auBerdem nur einmal eo 
berechnet werden. Der Aufwand zur Berechnung eines Koeffizienten fiber den Goertzel-Algorithmus ist in der 
praktischen Realisierung auf einem modemen Einchip-Signalprozessor TMS320C30 geringer als der ffir eine 
komplette FFT. (Fig. 17 zeigt den Rechenaufwand fur verschiedeneTransformationsiangen N). 

Fig. 18 zeigt dazu, daB innerhalb der Rechenzeit, die eine schneile Fourier-Transformation bendtigt, in vielen 
Fallen sogar einige Goertzel-Algorithmen gerechnet werden konnen. Wie aus beiden Fig. 1 7 und 18 hervorgeht, 55 
ist es beispielsweise bei einer Transformationslange von N = 512 gunstiger, bis zu 13 Spektralkoeffizienten fiber 
den Goertzel-Algorithmus einzein zu berechnen, als fiber eine schneile Fourier-Transformation. Bei der Anwen- 
dung des Goertzel-Algorithmus in der BARK-Transformation gewinnt man daher einen ganz betrachtlichen 
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GeschwindigkeitsvorteiL Der Rechenaufwand fur das in Fig. 12 angegebene Beispie! reduziert sich dadurch 
enorm; in den meisten Fallen braucht sogar nur ein Spektralkoeffizient berechnet zu werden. Die Berechnung 
einer schnellen Fourier-Transformation uber eine Koeffizientenanzahl, die nicht einer Potenz von 2 entspricht, 
ist zwar m6glich, fuhrt aber im allgemeinen zu einem unverhaltnismaBig hdheren Aufwand aufgrund der 
5 komplizierteren Struktur der dabei einzusetzenden Verfahren. 

Ein weiierer Vorteil des Goertzel-Algorithmus besteht demgegenQber darin, daB beliebige Frequenzrasterun- 
gen Qber dasselbe Verfahren (d. h. uber das gleiche Programm) gerechnet werden kdnnen. Dazu ist lediglich die 
LSnge des Eingangsvektors und der Winkel 2nk/N gemSB der zu berechnenden Linie k neu vorzugeben. Obwohl 
der Goertzel-Algorithmus als rekursives Filter am Stabilitatsrand arbeitet. treten dennoch keine Stabilitatspro- 
to bleme (z. B. Oberlauf des Rechenwerks) auf, da der Algorithmus nach der Berechnung eines jeden Koeffizienten 
der BARK-Transformation wieder neu aufgesetzt wird 

FQr die BARK-Transformation ergibt sich somit gegenuber dem in Fig. 12 gegebenen Beispiel ein weitaus (um 
den Faktor 224) geringerer Gesamtaufwand zur Ermittlung des gehorrichtigen Kurzzeitspektrums: 25 Goertzel- 
Algorithmen zur Berechnung der komplexen BARK-Koeffizienten der 25 Teilbander und 1 1 Polyphasenfenster 
15 zur Einstellung der 1 1 verschiedenen Bandbreitea 

Beide Teilalgorithmen erfordem im Kern ausschiieBlich Multiplikations-Additions-Operationen und sind 
daher fQr eine AusfQhrung auf einem Signalprozessor pradestiniert Die konsequente Anwendung des Prinzips 
der Reziprozitat und die Polyphasen-Filterung fuhrt fflr zu der in Fig. 20 fOr eine Abtastfrequenz von fa = 
44»l kHz angegebenen Betragsubertragungsfunktionen der Einzelfilter einer BARK-Transformation. aufgetra- 
20 gen Qber einer linearen Frequenzachse normiert auf die Abtastfrequenz von fa = 44.1 kHz. Die DARK-Trans- 
formation weist ein konstantes Obertragungsverhalten im DurchlaBbereich, eine ausgezeichnete Sperrdamp- 
fung und eine hervorragende Weitabselektion (keine Kebenmaxima) auf. 

In Fig. 19 ist der Frequenzgang der Summe der Amplitudenubertragungsfunktionen angegebender BARK- 
Filterkaskade angegeben. 

25 

Besonderheiten der BARK-Transformation 

Die BARK-Transformation besitzt in der hier vorgestellten Form viele ideale Eigenschaften. Die Nachteile 
herkdmmlicher Frequenzzerlegungen uber Transformations- und Teilbandverfahren kdnnen damit iiberwunden 
30 werden. Die Hauptvorteile der BARK-Transformation sind: 

— Gehdrrichtige Frequenzauflosung: Die BARK-Transformation folgt dem von Zwicker ermittelten Zu- 
sammenhang fiir die Frequenzgruppenbreite an der jeweiligen Mittenfrequenz. Der Fehler, der durch die 
Vorgabe eines Frequenzrasters bei der iterativen Ermittlung der Grenzfrequenzen der frequenzgruppen- 

35 breiten Teilbander entsteht, kann beliebig klein gehalten werden, wenn nur das Raster fein genug ist 

— Feinstmogliche Zeitaufldsung: Die Impulsaniwort eines (Band-)Filters kann ais Fensterfunktion fur eine 
Transformation aufgefaBt werden. Die Transformationskoeffizienten stellen im Frequenzbereich Reprasen- 
tanten fur die (Fourier-)Transformierten des Fensters im Zeitbereich dar. Deshalb muB sich bei der Varia- 
tion der "Bandbreite eines Koeffizienten" auch die Auflosung im Zeitbereich andera Bei der BARK-Trans- 

40 formation ist daher fur jede Frequenzgruppe das Produkt aus der Bandbreite B und der ZeitauflOsung T 
konstant und gleichzeitig minimal: B-T = 0^. Diese Grenze. die als Heisenberg-Grenze bezeichnet wird, 
kann von keiner anderen Frequenzanalyse unterschritien werden. Der Filterentwurf entlang der Heisen- 
berg-Hyperbel gewahrleistet, daB weder die Zeitauflosung noch die Bandbreite der Teilbander durch die 
BARK-Transformation verschlechtert wird. 

45 - Linearphasige Teilfilter: Aufgrund der Linearphasigkeit der Teilfilter und den gezeigten idealen Eigen- 

schaften der aus der Iteration resultierenden Betragsubertragungsfunktionen kann die BARK-Transforma- 
tion in dieser Form als bestens geeignet zur Nachbildung der Gehoreigenschaft angesehen werden, Fre- 
quenzgruppen zu bilden. Die linearphasigen Teilfilter erfordern keine Laufzeitentzerrung durch Allpasse. 
Da Qberdies die Filterlangen fur alle Teilfilter gleich sind, erQbrigen sich auch Laufzeitanpassungen durch 

50 VerzOgerungsglieder. 

— Rechteckffirmige Teilbander: Der iterative Filterentwurf Qber die Fenstermethode kann interaktiv an 
einem Rechner erfolgen. Der Entwurf gewahrieistet - trotz der Beschneidung der Spaltfunktion auf eine 
endliche Unge — Betragsubertragungsfunktionen fur die frequenzgruppenbreiten Teilbander, die sich 
durch einen konstanten DurchlaBbereich und eine gute Weitabselektion wegen der fehlenden Nebenmaxi- 

55 ma auszeichnea AuBerdem wird ein nahezu konstantes Obertragungsverhalten der gesamten Filterbank 
erreicht Zu den ausgezeichneten Eigenschaften, vor allem der hohen Flankensteilheit bei relativ geringem 
Realisierungsaufwand gegenfiber einzelnen Bandfiltem, tragt nicht zuletzt die Polyphasen-Filterung beL 

— Okonomische Realisierung: Bei der Realisierung wird ausgenutzt, daB die Frequenzgruppenbreiten im 
unteren Frequenzbereich nahezu konstant sind Es muB daher mit weit weniger Fensterfunktionen im 

60 Zeitbereich multipliziert werden, als Teilbander nach der Transformation entstehen. 

Die mit der BARK-Transformation gehindene Frequenzanalyse Ober nahezu rechteckfdrmige Teilbander 
bildet den idealen Ausgangspunkt fiir die weitere Entwicklung eines gehdrrichtigen MeBverfahrens einerseits 
und gehdrangepaBter Q)dierverfahren andererseits. Das Transformationsverfahren laBt sich dariiber hinaus 
65 auch auf andere Anwendungen Qbertragen. bei denen einen Transformation aquidistanter Abtastwerte vom 
Zeitbereich in den Spektralbereich benStigt wird, wobei die Spektralkoeffizienten auf der Frequenzachse 
nicht-aquidistant angeordnet sind und dennoch die Heisenberg-Bedingung erfullt bleiben soil. 
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Patentanspruche 

1. Frequenzanalyseverfahren 

bei dem aus einem zeitdiskret in Form von Abtastwerten vorliegenden Eingangssignal mittels Fensterfunk- 
tionen uberlappend aufeinanderfolgende BIdcke von Abtastwerten entnommen werden, worauf die so mit 
den Fensterfunktionen bewerteten. BIdcke zur Erzeugung eines Satzes von Koeffizienten furier-transfor- 
miert werden. 
dadurch gekennzeichnet, 

daB jeweils ein solcher Block des Eingangssignals mit einem Satz von unterschiedlichen. symmetrischen 
Fensterfunktionen gleicher L§nge bewertet wird, um einen Satz von mit den jeweiligen Fensterfunktionen 
gewichteten Bldcken von Abtastwerten zu erzeugen. wobei sich die Fensterfunktionen durch die unter- 
schiedlichen Bandbreiten ihrer Fourier-Transformierten voneinander unterscheiden, und 
daB jeder der zeitgleich erzeugten Blocke von gefensterten Abtastwerten eines Satzes einer eigenen 
Fourier-Transformation unterworfen wird» 

— indem zu jedem der zeitgleich erzeugten BIdcke wenigstens ein Koeffizient berechnet wird, dem die 
Bandbreite der Fourier-Transformierten dieser Fensterfunktion zugeordnet ist, und 

— indem die Koeffizienten so gewahit werdea daB die thnen zugeordneten Frequenzbinder sich, 
bezogen auf ihre Bandbreite, nur geringfugig uberlappen oder nur geringfugig voneinander beabstan- 
det sind. 

2. Frequenzanalyseverfahren nach Anspruch 1, dadurch gekennzeichnet. daB die Fensterfunktionen sinx/ 
x-Funktionen sind, deren Mulldurchgange in Abtastzeitpunkte von Abtastwerten gelegt sind, so daB sich 
eine Zeitauflosung entsprechend dem gegenseitigen Abstand der Mulldurchgange ergibt wobei die Zeitauf- 
Idsung umgekehrt proportional zur doppelten Bandbreite der Fourier-Transformierten ist 

3. Frequenzanaiyseverifahren nach Anspruch 2, dadurch gekennzeichnet, 

daS die Fourier-Transformation eine diskrete Schnelle Fourier-Transformation (FFT) ist, die eine der 
Anzahl von Eingangswerten im Zeitbereich entsprechende Anzahl von Koeffizienten im Frequenzbereich 
liefert. deren Abstand auf der Frequenzachse der Abtastfrequenz geteilt durch die Lange der Transforma- 
tion entspricht und deren Bandbreite durch die Fensterfunktion vorgegeben ist, und 

daB aus diesen Koeffizienten nur solche selektiert werden, die zusammen mit zeitparallel aus dem Satz der 
aus weiteren schnellen Fourier-Transformationen auf demselben Block von Abtastwerten mit unterschiedli- 
chen Fensterfunktionen zeitparallel zur Verfugung stehenden Koeffizienten bei anderen Mittenfrequenzen 
sich eine Kaskade von Bandfiltern aufbauen laBt, so daB die sich ergebenden Frequenzbander im wesentli- 
chen lilckenlos aneinander angrenzen. 

4. Frequenzanalyseverfahren nach Anspruch 2. dadurch gekennzeichnet, 

daB zur Verbesserung der Selektionseigenschaften der so entstehenden Kaskade von Bandfiltern die Lange 
der unterschiedlichen Fensterfunktionen eines Satzes auch ein Vielfaches der Transformationslange betra- 
gen kann. 

wobei vor der Transformation die Vielfachheit von Eingangswerten durch Aufteilung und additive Oberla- 
gerung (Superposition) der bereits gefensterten Eingangswerte wieder auf die Transformationslange redu- 
ziert wird(Polyphasenaddition). 

5. Frequenzanalyseverfahren nach Anspruch 2. dadurch gekennzeichnet, daB die selektierten Koeffizienten 
aus dem Eingangssignal mittels eines Goertzel-Algorithmus zweiter Ordnung berechnet werden. 

6. Frequenzanalyseverfahren nach einem der Anspruche 1 bis 5, gekennzeichnet durch 

ein iteratives Verfahren zur Bestimmung der diskreten Mittenfrequenzen und Bandgrenzen zur Festlegung 
der Fensterfunktionen, welches folgende Schritte umfaBt: 

— Setzen einer Lauf frequenzvariablen (fmiiie) auf eine Startfrequenz (foben)* die auf dem durch die 
Lange der Fourier-Transformation gewahlten Frequenzraster liegt; 

— Verschieben der Lauffrequenzvariablen (fmitte) auf dem Frequenzraster um einen Schritt (df); 

— Ermittlung der sich aufgrund einer gegebenen, nicht linearen Funktion der Mittenfrequenz ergeben- 
den Bandbreite; und 

— Wiederholung der Iteration solange, bis eine Bandgrenze des zu der Lauffrequenzvariablen (fmitte) 
gehorenden Bandes mit einer Bandgrenze eines vorher ermittelten Bandes zusammenfallL 

7. Frequenzanalyseverfahren nach Anspruch 6, dadurch gekennzeichnet, daB die gegebene, nichtlineare 
Funktion der Zusammenhang zwischen der Mittenfrequenz und der Frequenzgruppenbreite ist 
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Vorgaben: 
Abtastfrequenz: fa 
Frequenzraster: #lines 



J 




foben = fa/2 
df = foben / #lines 



fmitte = foben 



L 



fmitte = fmitte - df 




funten = fmitte - (foben - fmitte) 



T 



BW = (foben - funten) 



foben = funten 




Ausgabe an Filterentwurf: 
BW(fmitte) 
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Audio signal frequency analysis method - using window functions to provide sample signal 
blocks subjected to Fourier analysis to obtain respective coefficients. 
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The frequency analysis method involves dividing an input signal, represented by discrete sample 
values, into overlapping blocks using a window function and subjecting them to a Fourier 
transformation, to obtain a set of coefficients. Each input signal block is evaluated via a set of 
differing symmetrical window functions of equal length, their Fourier transformations having 
different bandwidths. 

The simultaneously generated blocks are used to provide a set of characteristic Fourier 
transformations, with a coefficient provided for each block, dependent on the window function 
Fourier transformation bandwidth. The selected coefficients have frequency bandswhich are only 
slightly overlapping, or only slightly spaced apart. 

USE/ AD VANTAGE - Accurate simulation of human hearing characteristics. At least first five 
spectral coefficients processed together, e.g. in fast fourier transformation. 

Dwg.6/20 

US 5583784 A 

A frequency analysis method that simulates perception characteristics of human hearing, comprising 
the steps of: 

extracting from an input signal a sequence of blocks of discrete sampled values, wherein successive 
blocks of the sequence overlap each other in time; 

concurrently windowing one of the successive blocks using a plurality of differing, 
symmetricalwindow functions of substantially identical length, each window function having a 
corresponding Fourier transform with a different bandwidth, to produce a plurality of windowed 
blocks of weighted sampled values; and 

performing a Fourier transformation on each windowed block wherein for each windowed block at 
least one coefficient is generated having assigned thereto the bandwidth of the Fourier transform of 
the window function used to produce the windowed block, and wherein the at least one coefficient for 
each windowed block is chosen such that bandwidths assigned to the chosen coefficientsof the 
plurality of windowed blocks adjoin one another approximately without gaps and in a manner 
approximately free from overlapping. 

(Dwg. 12/20 
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